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ABSTRACT 

We report on the detection by Swift of GRB 080913, and subsequent optical/near-infrared follow-up 
observations by GROND which led to the discovery of its optical/NIR afterglow and the recognition 
of its high-z nature via the detection of a spectral break between the i' and z' bands. Spectroscopy 
obtained at the ESO-VLT revealed a continuum extending down to A = 9400 A, and zero flux for 
75001 < A < 9400 A, which we interpret as the onset of a Gunn-Peterson trough at z=6.695±0.025 
(95.5% conf. level), making GRB 080913 the highest redshift GRB to date, and more distant than the 
highest-redshift QSO. We note that many redshift indicators which are based on promptly available 
burst or afterglow properties have failed for GRB 080913. We report on our follow-up campaign and 
compare the properties of GRB 080913 with bursts at lower redshift. In particular, since the afterglow 
of this burst is fainter than typical for GRBs, we show that 2 m-class telescopes can identify most 
high-redshift GRBs. 

Subject headings: gamma-rays: bursts - radiation mechanisms: non-thermal - early Universe 
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1. INTRODUCTION 

The potential of gamma-ray bursts (GRBs) as bea- 
cons to the distant universe has long been recognized. 
The immense luminosity of both the prompt gamma- 
ray emission and X-ray and optical afterglows indicates 
that GRBs should, with present-day technology s uch as 
Swift /BAT be visible out to distanc es of z > 10 ( Lamb 



the 



& Reichart 112000 Gou et al. 20041. Due to their con 



nection to the death of massi ve stars (|Woosley 
Galama et al. |[1998 Hjorth et al. 
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lo ng GRBs probe the evolution of cosmic star 
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re-iomzation of the in 
tergalactic medium (Miralda-Escude 1998, Totani et al. 
2006), and the metal enrichment history of the Universe 



e.g. [Hartmann et al. |2004j Fynbo et al. 



2006; Berger et al. |2006{ |Prochaska et af 



2006 



2007) 



Savaglio 



Despite 

this potential and the high detection rate now delivered 
by Swift, finding bursts out to the highest redshifts has 
proved a challenging task. The hitherto highest redshift 
burst, GRB 050904 at z = 6.29 flKawai et al. |2006||Hais- 



lip et al 

that time oriTy~GRB 060927 at z 



20061, held the record for th ree years and in 
5.47 ( Ruiz- Velasco et 



al^[2007| came close. 



n part this owes much to the intrinsic difficulty in lo- 
cating such bursts. At z > 5.5 afterglows become essen- 
tially invisible to observers in the i?-band, where much 
follow-up is attempted. Often the mere time required 
to obtain a photometric selection of high-z candidates 
is so long that the afterglow is too faint for spectro- 
scopic follow-up. Furthermore, although the campaigns 
on GRB 050904 were extremely successful, this burst was 
far from a typical event. Indeed, its peak optical lu- 
minosity rivalled those of the exceptional bright GRBs 
990123 and 080319B (Akerlof et al. 1999; Kann et al. 
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Not corrected for Galactic foreground reddening of E(B-V) = 0.043. Converted to the AB system for consistency with Fig. I2I using 
the following coefficients for the J band: HAWK-I: +0.98 mag, Gemini-N: +0.96 mag, ISAAC: +0.96 mag, GROND: +0.91 mag, SOFI: 
+0.96 mag, Subaru: +0.94 mag. 



2007a, Racusin et al. 2008). To identify more "typi- 
cal" afterglows at high redshift is a much more arduous 
task, requiring rapid response, multi-color observations 
and ultimately rapid and deep spectroscopy. Such situa- 
tions must inevitably be somewhat fortuitous, requiring 
a burst to be visible to large telescopes almost immedi- 
ately after its occurrence, with good weather and appro- 
priate instrumentation. This perhaps explains, at least 
in part, why bursts more distant than GRB 050904 have 
been extremely difficult to find. Yet, the present detec- 
tion rate of GRBs at z ^5 is about what is predicted on 
the basis of the star formation rate, so the lack of many 
high-z GRBs may well be due to the fact that there are 
few of them. 

Here we report the d iscovery of GRB 080913 with Swift 



(Gchrels et al. 2004), and subsequent follow-up obser- 



vations which identify this GRB to have originated at a 



redshift of z — 6.7, the highest known to date. In Sect. [2] 
we describe the observational effort in all wavelengths 
from hard X-rays to near-infrared, and present the re- 
sults of these observations. In Sect. [3] we discuss several 
aspects of our observational findings. 

2. OBSERVATIONS AND ANALYSIS RESULTS 

2.1. Swift BAT, XRT and UVOT measurements 

Swift/BAT triggered on GRB 080913 (t rigger 324561) 
on Sep. 13, 2008 at T n = 06:4 6:54 UT (|Schady et al. 



on scp 
2008 



^ Stamatikos et al. ||2008| . The BAT light curve 
shows multiple overlapping peaks with a Tg duration 
(the time interval during which 90% of the fluence is 
measured) of 8±1 s. The peak count rate was 800 
counts/s (15-350 keV). Using a 64 ms binned light cur ve, 



we compute spectral lags (e.g. |Hakkila et al. 2007) of 
-lT 4+0.124 s for the 100-150 keV vs. 50-100 keV band, 



GRB 080913 
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and Q.148±g;g|| s for the 100-150 keV vs. 15-50 keV 
band, consistent with an independent estimate by Xu 
(2008). 

The time-averaged spectrum from To — 3.8 s to To+5.2 
s can be adequately fit by a power law with an exponen- 
tial cutoff. This fit gives a photon index of 0.46±0.70, 
and E peQfc = 93 ± 56 keV (x 2 / dof = 38.5/56). The total 
fluence in the 15-150 keV band is (5.6 ± 0.6) x 10" 7 
erg/cm 2 , and the 1-s peak flux measured at To+0.11 s 
in the 15-150 keV band is 1.4 ± 0.2 ph/cm 2 /s. A fit to 
a simple power law gives a photon index of 1.36 ± 0.15 
(x 2 /dof = 44.6/57) (all the quoted errors are at the 90% 
confidence level) . A combined fit of the Swift /B AT (15- 
150 key) and Konus-Wind (20-1300 keV) data (pal'shin 



GROND, a simultaneous 7-channel imager ( |Greiner"et 



al. 2008a I mounted at the 2.2 m MPI/ESO telescope at 



et al. 20081 in the time interval T - 4.1 to T +4.7 s 



using a Band function (Band et al. 1993; /3=-2.5 fixed) 
yields a = -0.82^°;^ and E peak = 121±§| 2 keV (where 
Epeafc is the energy at which most of the power is emit- 
ted, and a and (3 are the low- and high-energy photon 
indices, respectively), and an energy fluence in the 15- 
1000 keV band for th e 8.8 s interval of 9xl0 -7 erg/cm 2 



(iPal'shin et al. 1120081) 

" -II. I I ' 



Swift slewed immediately to the burst and the X-ray 
Telescope (XRT, Burrows et al. 2005a) began its au- 
tomated observing sequence at 06:48:30UT, 96 s after 
the trigger. A fading X-ray afterglow was detected; the 
UVOT-enhanced X-ray position (Beardmore et al. 2008) 
is RA (J2000.0) = 04 h 22 m 54?66, Decl. (J2000.0) = - 
25°07'46 / .'2, with an uncertainty of l'/9 (radius, 90% con- 
fidence; see Fig. [TJ. 

The X-ray light curve (Fig. [2]) shows an initial decay 
rate (F cx t~ a ) of a - 1.2 from -100 s to -300s and 
from —400 to —1100 s, with a likely small flare with 
the same decay slope afterwards. At around To+1.8 ks 
(AT/T — 0.3), there is a substantial flare with a factor 
of — 5 increase in count rate. The evolution thereafter 
is difficult to characterize due to the sparse coverage: it 
could be a continued decay from 2 ks to —100 ks at the 
same slope of 1.2 but with an offset suggestive of energy 
injection, followed by a plateau, or it can be described as 
decaying at a slope of about 1.2 (Beardmore & Schady 
2008) from 1200 s to 10 6 s with two flares super-imposed. 

The X-ray spectrum (using the Swift software pack- 
age v2.9 and its associated set of calibration files) from 
To + 108 s to 1920 s is well fit by an absorbed power 
law of photon index T = 1.66 ± 0.14. A separate fit of 
the early flare and preflare times gives the same pho- 
ton index (1.69±0.25). The column density is consis- 
tent with the Galactic value of 3.2 x 10 20 cm -2 in the 
direction of the burst (Kalberla et al. 2005). The ob- 
served 0.3 — 10.0 keV flux at this time was 8.0^' o x 
10~ 12 erg cm -2 s _1 . This corresponds to an unabsorbed 
flux of 8.5tJ;J X 10- 12 ergcm- 2 s- 1 (0.3 - 10.0 keV). 

UVOT observations in white-light (100 s exposure) 
started 105 s after the BAT trigger, and subsequently all 
UVOT filters were used. The afterglow was not detected 
in any of the UVOT filters and the 3 a limiting white 
magnitude for the first finding chart exposure is >20.92; 
see Oates & Schady (2008) for the detailed upper limits 
in each filter. 

2.2. Optical/NIR photometry 
2.2.1. Observations 



La Silla (Chile), started observing the field at 06:52:57 
UT, about 6 min after the GRB. The imaging sequence 
began with 66 s integrations in the g'r'i' z' channels with 
gaps of about 16-18 s due to detector read-out and preset 
to a new telescope dither position. After about 15 min, 
the exposure time was increased to 115 s, and after an- 
other 22 min to 375 s. In parallel, the three near-infrared 
channels JHKs were operated with 10 s integrations, 
separated by 5 s due to read-out, data-transfer and K$- 
band mirror dithering. At 09:46 UT, the start of nautical 
twilight (Sun is 12° below horizon), GROND switched to 
the "NIR-only" mode, in which the CCDs of the 4 visual 
channels are switched off, and only imaging in JHKs is 
performed. Observations finally stopped completely at 
10:07 UT. Further GROND imaging was performed on 
Sep. 15, 07:12-09:28 UT, and Sep. 16, 06:38-09:47 UT. 

Imaging was also secured immediately after the burst 
and on the subsequent nights at the VLT (ESO) with 
the Focal Reducer and Spectrograph FORS in filters 
BVRIz, and with ISAAC and HAWK-I in J. Further 
imaging was obtained with NTT/SOFI, Subaru/IRCS, 
and Gemini-N/NIRI in the J band (Tab. [TJ. 

Data reduction was done in a standard way using 
IRAF routines. Photometric calibration of the GROND 
g',r',i',z' bands was performed using the Sloan spec- 
trophotometric standard star SA95-142, which was ob- 
served shortly after the GRB 080913 field. Magnitudes 
are therefore given in AB magnitudes since this is the 
natural photometric system for GROND (see Greiner et 
al. 2008a). In order to match the different z'J-band 
filters used, the field calibration of the GROND filter 
magnitudes was compared against the VLT (or Gem- 
ini/Subaru) magnitudes using —50 field stars. For in- 
stance, the rms scatter between the GROND z' and 
FORS2 Gunn-z was < 0.06 mag, so no color transfor- 
mations were applied. Calibration of the field in JHK 
was performed using 2MASS stars which were chosen as 
close as possible to the GRB afterglow (and with a nicely 
sampled point-spread function), to reduce the error in 
the calibration due to flatfielding, which is (for GROND) 
0.5-1.0% on small scales and around 2% over the whole 
array. The magnitudes of the selected 2MASS stars were 
then transformed into the GROND filter system and fi- 
nally into AB magnitudes using J(AB) = J(Vega) + 
0.91, -ff(AB) = ff(Vega) + 1.38, K(AB) = ^(Vega) + 
1.81 (for details, see Greiner et al. 2008a). In this way a 
set of nine secondary standard stars was created (Tab. [21. 

2.2.2. The photometric redshift estimate 

After stacking the GROND exposures taken during 
the first 16 min (06:53-07:07 UT) the G ROND pipeline 



reduction (Kiipcii Yoldas, et al. 2008) found a faint 



source at RA (J2000.0) = 04 n 22 m 54!74 Decl. (J2000.0) 
= -25°07'46'.'2 (0'.'5 error), which was only detected in 
the z' JHKs-b&nds, but not in shorter- wavelength bands 
(Figs. JTJJ3|, and thus suggested a redshift above 6 (Rossi 
et al. 2008; Greiner et al. 2008b). This was confirmed by 
VLT/FORS2 RzIVBR imaging, starting 45 minutes af- 
ter the burst, the only detection bein g in the z band with 
z(AB) = 23.1 ( |Vreeswijk et al 120 08 ) . A fit to the simul- 
taneously obtained 7-filter GROIN D spectral energy dis- 
tribution (SED), using Hyper-z (Bolzonella et al. 2000) 
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TABLE 2 

Local photometric standards; the JHKs magnitudes are in the Vega system, the g'r'i'z' magnitudes in AB. 



No 


Coordinates (J2000) 


g' 


r' 


•/ 

i 


z 


J 


H 


K s 


1 


04 h 22 m 56!l 


-25°07'10" 


_ 


20.52±0.02 


19.284=0.01 


18.724=0.01 


17.504=0.04 


16.90±0.03 


16.90±0.06 


2 


04 h 22 m 54.4 


-25°07'27" 










20.264=0.05 


19.314=0.06 


18.61±0.13 


3 


04 h 22 m 53=9 


-25°07'37" 










19.97±0.05 


18.914=0.05 


18.104=0.11 


4 


04 h 22 m 55.2 


-25°07'39" 


16.974=0.01 


16.14±0.01 


15.82±0.01 


15.654=0.01 


14.76±0.02 


14.25±0.03 


14.204=0.04 


5 


04 h 22 m 55.3 


-25°07'49" 










20.01±0.05 


18.86±0.06 


18.574=0.14 


6 


04 h 22 m 54.2 


-25°07'57" 


15.964=0.01 


15.57±0.01 


15.39±0.01 


15.29±0.01 


14.54±0.02 


14.18±0.03 


14.194=0.04 


7 


04 h 22 m 55=8 


-25°08'04" 


18.73±0.01 


17.26±0.01 


16.22±0.01 


15.76±0.01 


14.57±0.02 


13.88±0.03 


13.70±0.05 


8 


04 h 22 m 57!4 


-25°08'11" 




20.904=0.04 


20.20±0.03 


19.81±0.04 


18.34±0.05 


17.43±0.04 


16.714=0.07 


9 


04 h 22 m 52.7 


-25°08'36" 




20.044=0.02 


18.71±0.01 


18.15±0.01 


16.904=0.04 


16.24±0.03 


16.174=0.05 



results in a photometric redshift of z = 6.44±0.30 (Fig. 
[3| Greiner et al. 2008b). 

2.2.3. The afterglow light curve 

The optical/NIR light curve (Fig. [2| during the first 
few hours is described by a power law of 0.98±0.05. The 
third GROND data point at 2000 s post-burst is simulta- 
neous to the X-ray flare, and is clearly above that power 
law decay, suggesting that we see enhanced optical emis- 
sion from that X-ray flare (Fig. [2]). If we ignore this third 
data point, the reduced \ 2 improves substantially, and 
the new power law decay is 1.03±0.02. After the GRB 
location became visible again after 2 days (clouds pre- 
vented observations on Sep. 14), at ~180 ks post burst, 
the optical/NIR brightness as measured by GROND is 
nearly identical to that at 10 ks, suggesting a prolonged 
plateau phase (Fig. [2]). Later monitoring with GROND, 
VLT, NTT, Gemini^N and Subaru shows pronounced 
variability which can be, in principle, fit by three log- 
Gaussian flares. At the same time, also the X-ray light 
curve shows an enhanced flux level. It is likely, though 
it cannot be proven, that the X-ray light curve also con- 
sisted of flares, which are not resolved due to the faint- 
ness and the somewhat poorer sampling as compared to 
the optical/NIR band. In any case, the enhanced X-ray 
flux lets us prefer that the optical lightcurve is probably 
a flare rather than a plateau. Our J band data point 
(J = 24.61 ± 0.13) at 8 x 10 5 s is again well above a 
power law extrapolation after the two flares. 

2.2.4. The GROND- derived afterglow broad-band spectrum 

We fit the NIR SED (the z' and Gunn-z bands are 
affected by Ly-a and are not included) with no extinction 
and with three different dust models (LMC, MW, SMC). 
Here we assume that the afterglow spectrum should have 
a power law shape, with F cx (note that the spectral 
index /3 is related to the photon index T from the X-ray 
spectral fitting by (3 — T - 1). The best-fit power law has 
j3 = 1.124=0.16. No evidence for extinction in the host 
frame is found. 

2.3. Optical spectroscopy 

Immediately after determining the photo-z, we trig- 
gered our target-of-opportunity program at ESO (PI: J. 
Greiner) , and obtained an optical spectrum of the after- 
glow of GRB 080913 in the 7500-10500 A (grism 600z) 
region with FORS2/VLT on Sep. 13, 2008. Only one 



1800 s and one 600 s exposures were possible before dawn. 
The spectrum revealed a continuum disappearing blue- 
wards of a break around 9400 A. Interpreting this break 
as the onset of the Lyman-a fore st a preliminary re dshift 
estimate of z — 6.7 was inferred ( Fynbo et al. |2008b[ ). A 
second set of exposures (7x1800 s) was acquired on Sep 
16, 2008, with identical settings. The first night's spec- 
trum is shown in Fig. |4j The afterglow is also detected 
in the spectrum from 2 nights later, but with smaller 
signal-to-noise ratio. 

In order to constrain the redshift we first tried to 
search the spectrum for a significant metal absorption 
line. Unfortunately, no such line seems to be present 
in the spectrum. Therefore, the redshift can be based 
on the shape of the break at 9400 A only: the flux is 
zero in the interval 7500-9400 A, and clearly non-zero 
above 9400 A, strongly supporting the detection of the 
Lyman break, i.e. a high-z object. We first tried a simple 
cross-correlation analysis using the "fxcor" task in IRAF. 
As templates we used the spectra of GRBs 050730 (z = 
3.969) and 060206 (z = 4.048, redwards of the center 
of the DLAs of these spectra) . This analysis leads to a 
redshift of z = 6.69 4= 0.02. 

In order to make a more refined analysis we follow 
Totani et al. ( 2006 1 and fit the spectrum using two dif- 
ferent assumptions about the origin of the break. The 
damping wing can either be interpreted as that of a 
Damped Lyman-q (DLA ) absorber in the host galaxy 
( Jakobsson et al. 2006b) or that of IGM neutral hydro- 
gen ( |lVladau &; Rees ||2000[ ). We performed a chi-square 
analysis with the three parameters of N^a (column den- 
sity of DLA), xm (neutral fraction of IGM hydrogen), 
and z (redshift of the GRB). We assume that the DLA 
has the same redshift z, and the neutral IGM hydrogen 
is distributed from Zigm,i to z, meaning that there is 
no ionized bubble around the host galaxy. We assume 
z igm,i — 6.0, but this has little influence on the fit if 
z — Zigm,1 >~ 0-3 (Totani et al. 2006). No prior is 
assumed for these model parameters. We obtained the 
best-fit at (logiVHi/cm- 2 ,:ziHi,z) = (19.74,1.0,6.71), in 
which the damping wing is dominated by IGM. On the 
other hand, a DLA dominating case of (21.41, 0.001, 
6.68) is also consistent with the data (1.2 a deviation 
from the best fit). Marginalizing TVjji and xjji, we ob- 
tain the 95.4% confidence limits on z as 6.67 < z < 6.72, 
where the DLA and IGM are dominant in the lower and 
upper limits, respectively. 

When we fix xm — 0.001 (almost fully ionized IGM 
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with negligible effect on the damping wing), we obtain 
95.4% limits on N m as 20.29 < logA^i/cm" 2 < 21.41 
marginalizing the redshift, with the best fit value of 
log -/Vm/cm -2 = 20.99 at z = 6.69. When we assume 
that TVhi is negligibly small, we obtain a 95.4% lower- 
limit of xhi > 0.35 again marginalizing redshift. This 
means that, if the column density of the DLA in the host 
galaxy is sufficiently small to have a negligible effect on 
the observed damping wing, we need a significantly neu- 
tral IGM indicating that the re-ionization has not yet 
been completed at z — 6.7. Although this possibility 
is completely dependent on the assumption about DLA 
column density, it is not in contradiction with the 95% 
upper bound of xm < 0.60 at z = 6.3 obtained from GRB 
050904 (Totani et al. 2006). A significantly high neutral 
fraction of IGM at z > 6.7 has also been implied from 
evidence for the GP damping wing in the spectra of z=6 
quasars (Mesinger & Haiman 2004, 2007) and the lumi- 
nosity function evolution of Ly-a emitters (Kobayashi et 
al. 2007; Ota et al. 2008). These claims remain contro- 
versial, since different sightlines to quasars can produce 
large variations in the Ly-a absorption spectra (Bolton 
& Haehnelt 2007) , and the apparent evolution of the Ly- 
a luminosity function may be explained by the evolution 
of the mean IGM density alone (Dijkstra et al. 2007). 

2.4. The X-ray to optical SED 

The combined XRT and GROND/VLT spectral energy 
distribution was fit (except the z' band since it is affected 
by Ly-a) in three separate time intervals: the early de- 
cay except the second flare at To+2000 s, i.e. the interval 
360-1200 s, the intermediate brightness level (5.16-12.6 
ks), and at very late times (16.7-550 ks). We used a 
single as well as a broken power law, modified by dust 
extinction in the rest-frame FUV (zdust in XSPEC) for 
Milky Way, SMC and LMC, and hydrogen absorption 
at X-rays. The hydrogen absorption was fit for all the 
data, and then fixed at that value for the individual in- 
tervals. For the single power law fit, we obtain a photon 
index of 1.7±0.1 and a marginal reddening E(B-V) of 
0.096±0.045 (x 2 /dof = 8.4/10) assuming LMC dust in 
the GRB host (see Fig. M. The results for SMC or 
MW dust are the same within the errors. Forcing the 
dust extinction to be zero results in a somewhat worse 
fit (x 2 /dof = 15/11). We note that the GROND-XRT 
SED for this first time interval likely is affected by en- 
hanced X-ray emission from the early flare (at 300-700 
s), thus resulting in a flatter spectral slope than derived 
from the GROND data alone. 

Using a broken power law fit instead and fixing the 
XRT power law slope to 1.7 (see above), we obtain a 
photon index of 2.1±q[qI for the GROND data (corre- 
sponding to the energy index of 1.1; see above) and neg- 



ligible dust, i.e. E(B-V) = 0.026 



+0.027 
-0.054- 



This implies that 



the early X-ray afterglow emission is contaminated by 
the flare emission, leading to a flatter X-ray-to-optical 
SED. The goodness of the fit is not better than the sin- 
gle power law fit (% 2 /dof = 11/9), so a broken power law 
is anyway not required by the data. For the two late- 
time SEDs, the uncertainties are larger due to the lower 
signal-to-noise ratio, and thus are equally well fit by a 
single power law. 

2.5. Late-time XMM-Newton observation 



A triggered XMM-Newton observation was performed 
at 380 ks post-burst. The afterglow is clearly detected. 
There is no evidence for any variability in the XMM 
lightcurve. The spectrum is acceptably fit with a power- 
law with fixed Galactic absorption (3.2x 10 20 cm -2 ). The 
best-fit photon index is V = 2.0±0.2. The unabsorbed 
flux is 2.2xl0 -14 erg/cm 2 /s in the 0.3-10.0 keV band. 

3. DISCUSSION 

With z=6.695±0.025, GRB 080913 is the most distant 
burst detected so far (and actually the second most dis- 
tant spectroscopically confirmed object after a galaxy at 
z = 6.96, as well as the most distant X-ray sour ce) for 



which a redshift has been determined (lye 20081. In a 
concordance cosmology model (Spergel et al. 2003) , with 
H = 71 km/s/Mpc, Q M = 0.27, A = 0.73, the corre- 
sponding luminosity distance is 67 Gpc (2 x 10 29 cm), 
and the age of the Universe at that time is 825 Myr (6% 
of the present age) . 

3.1. Comparison to other GRBs 
3.1.1. GRB properties 

With the above numbers, the isotropic energ y release 
is E ?:so » 7 x 10 52 erg (1 keV - 10 MeV; see also |Pal'shin" 
et al. |20"08| . This is typical of the long GRB population 
at large. Compared to the properties of GRB 050904, 
the other previous GRB at z>6, GRB 080913 has a sub- 
stantially shorter duration, lower gamma-ray luminosity 
as well as much dimmer (~5 mag at early times) after- 
glow. Thus, it appears that gamma-ray bursts at this 
early epoch show the same large diversity as the low-z 
bursts. 

3.1.2. On the burst duration 

The observed bimodal hardness-duration distribution 
of the prompt emission from GRBs has led to the dis- 
tinction J3e_tweejij£n|£ (<^s)a^ 



bursts (|Mazets et al. [T 
liotou etal ||1993| ). L 



1981 Norris et al. 1984 Kouve- 



ong-duration gamma-ray bursts 
are thought to arise in jets created by the collapse of a 
massive star, short-duration bursts have b een suggested 
to emerge from a compact binary merge r (Eichler et al 



1989 Paczyhski 1987 Woosley 19931. GRB 080913 



is, m its rest frame, relatively short and hard, with a 
rest frame duration around 1 s, and rest frame E pea £; 
~ 990 keV for a Band function fit (Pal'shin et al. 2008). 
Other moderately high redshift bursts with rest-frame 
duration shorter than 2 s are GRBs 060206 (T 90 = 7±2 
s, z = 4.0), 051016B (T 90 = 4±0.1 s, z = 0.94) , 050406 
(T 90 = 5 ± 1 s, z = 2.44), 050416A (T 90 = 2.4 ± 0.2 
s, z = 0.6535), 000301C (T 90 ~ 2 s, z = 2.03), 040924 
(T 90 - 1.5 s, z = 0.859) and 050922C (T 90 = 5±1 s, 
z = 2.2) (Levan et al. 2007), as well as GRBs 060223: 
(T 90 = 11±2 s, z = 4.41), 060926: (T 90 = 8.0±0.1 s, z 
= 3.20), 070506: (T 90 = 4.3±0.3 s, z = 2.31), 071020: 
(T 90 = 4.2±0.2 s, z = 2.145), 080520: (T 90 = 2.8±0.7 
s, z = 1.545). These bursts all show several of the typ- 
ical signs of a massive stellar origin, e.g., star-forming 
host galaxies, a gas-rich line-of-sight, and (in the case of 
GRB 050416A) an associated supernova. Since the long- 
to short-duration separation at ~2 s was derived in the 
observer's frame for BATSE bursts for which redshifts 
are not available in general, we consider it most plausi- 
ble that GRB 080913, despite its intrinsically relatively 
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short duration, is a member of the long-duration popu- 
lation of GRBs. Also, GRB 080913 is compatible with 
the previously known lag-luminosity correlation of long- 
duration bursts (, and is also consistent with the Amati 
relation. While a merger origin for GRB 080913 can not 
be definitely ruled out, the short duration also leads to 
an interesting question for the collapsar scenario: How 
can a massive star produce a burst as short as 1 second? 

3.1.3. Afterglow properties 

In order to compare the afterglow of GRB 080913 with 
those of previous GRBs we analyze its intrinsic properties 
both in the optical as well as the X-ray regime within 
large afterglow samples obtained for other GRBs. Using 
= 1.1 we sh ift the light curve to z = 1 following |Kann| 



et al. | (J2006J) . The result is presented in Fig. |6| and 
compared with other GRB afterglow light curves. At 
early times, the afterglow is fainter than the mean of 
the sample, and especially much fainter than the optical 
afterglows of the GRBs with the second and third highest 
redshifts, GRB 050904 and GRB 060927, respectively. 
At one day after the GRB at z — 1, following a strong 
rcbrightcning, R c = 18.80 ±0.13 and M B = -24.1 ±0.2, 



which i s brigh ter than the mean of the sample of Kann 
et al. | ( |2Q07b] ), = -23.0 ± 0.4, but not exceedingly 
so. Indeed, even after the strong rebrightening, there are 
several afterglows that are considerably brighter at this 
time in the same frame. 

Fig. [7] shows the luminosity evolution of the X-ray 
afterglow of GRB 080913 in comparison to the popula- 
tion of 110 X-ray afterglows of Swift bursts with known 
redshift detected by July 2008. The original data were 
obtained from the X-ray light curve repository (Evans et 
al. 2007) . Luminosities were calculated following Nousek 
et al. (2006). The spectral slope (3 of the SED was 
obtained by fitting an absorbed power-law, F(E \> S ) ~ 
exp(-A^ al a(E ohs )) exp(-iV£ ost a(E host )) u~P with 
-Ehost = £obs (1 + z), while fixing the Galactic column 
density to the value given by Kalberla et al. (2005). 
In Fig. [7] the low luminosity end is mainly occupied 
by short burst afterglows, which on average separate 
clearly from their long burst cousins. Some short bursts 
are indicated (GRBs 051221A, 060313, 060121 at z=4.6: 
de Ugarte Postigo et al. 2006, or GRB 070714B) as 
well as the second and third most distant bursts GRB 
050904 (Haislip et al. 2006; Kawai et al. 2006) and 
GRB 060927 (Ruiz-Velasco et al. 2007). These events, 
together with GRB 080913, indicate a rather broad 
luminosity distribution even at high redshift. Fig. [7] 
makes clear that similar to the optical/NIR bands, the 
rest-frame X-ray afterglow of GRB 080913 is rather 
faint. Compared to the long burst ensemble it belongs 
to the low luminosity members, while compared to the 
short burst ensemble it occupies the high-luminosity 
region. 

3.2. Modelling of the afterglow 

The temporal and spectral power law indices of the 
early (t < i bl ~ 10 4 s) optical/NIR afterglow, a optj i = 
1.03 ± 0.02 and (3 op t = 1.12 ± 0.16, are consistent with 
the closure relation a pt.i — l-5/3 pt ~ —0.5, which cor- 
responds to the standard model with the optical band 
above the cooling frequency v c and the typical syn- 
chrotron frequency v m . The electron energy distribution 



index p ~ 2.2, inferred from p = (4a op t,i + 2)/3 and/or 
p = 2/3 opt of the standard model prediction (Meszaros & 
Rees 1997; Sari et al. 1998), is identical to the canonical 
one. 

The early optical/NIR data can also constrain the val- 
ues of other parameters (e.g., initial isotropic kinetic 
energy Ei, energy fraction in electrons e e and that in 
magnetic fields es) of the forward shock model. If the 
circum-burst environment is interstellar medium, the re- 
quirement of v op t > max(y m , j/ c ) at t = 560 s and the 
H band flux density F VH ~ 4 /xJy at t = tf,i, can 
be translated to a set of loose constraints. In general, 
10~ 5 < 6b < e e < 1, £^,53 < 1.0 and/or a density of 
n > 100 cm -3 can fit the afterglow, which is quite reason- 
able except maybe for the relatively high density. If the 
circum-burst environment is a free wind, then the same 
observational constraints would require a more ad hoc 
set of parameters, i.e. 10~ 5 < es < e e < 1, £^53 < 1.0 
and A* > 1.0 (where A* = A/(5 x 10 u gcm _1 ), and A 
is the normalization in the density profile p = Ar~ 2 ). 
While a high-density medium is expected at high-z (e.g. 
Gou et al. 2007), a large wind parameter A* > 1.0 is not 
expected in view of the low stellar metallicity at high-z. 
We therefore conclude that the afterglow modeling fa- 
vors a constant density medium, although a stellar wind 
medium is not ruled out. 

After the early decay, the afterglow intensity enters a 
plateau/flares phase. Due to the lower signal-to-noise 
ratio the optical to X-ray SED with a photon index of 
1.94±0.20 now has larger error bars, but is still consistent 
with Pox ~ 1.0 which indicates that the emission of both 
bands possibly has the same origin. 

For the flares interpretation (see solid line in Fig. j2|, 
the late time optical/NIR light curve is fit with log- 
Gaussian functions superimposed to the initial unbro- 
ken power law decay with temporal index fixed to be 
c*opt,i = 1-03. The best-fit amplitude and mid-time 
of the optical flares are F Ujp = 2.9 ±1.2 fj,3y, t p = 
(1.33±0.17) x 10 5 s for the first flare, F VjtP = (0.95±0.28) 
^Jy, t p — (3.00 ± 0.43) x 10 5 s for the second flare, and 
Fuj, P = (0.52 ± 0.05) /iJy, t p = (7.7 ± 0.4) x 10 5 s for 
the third flare. The late time X-ray flare happens al- 
most simultaneously with the optical flares. These late 
flares together with the early two X-ray flares (at To + 
300 s and To + 2000 s) are very likely the emission from 
late internal shocks, as already detected in many previ- 
ous Swift GRBs, long or short, and in GRB 050904 at 
comparable redshift (Burrows et al. 2005b; Barthclmy 
et al. 2005; Watson et al. 2006; Cusumano et al. 2006). 
The central engine of GRB 080913 in this case has accel- 
erated and emitted relativistic ejecta intermittently for 
a whole period of at least ~4x 10 4 s in its rest frame. 

For the plateau interpretation (see dotted lines of Fig. 
2), the late time optical light curve can be fit by a broken 



power law with a op t,: 



-0.19Io;ig during the plateau 



0.09 



until t = t b2 = 1.16t£$ x 10 5 s, and a 3 = 0.92+^ 
after the plateau. A joint fit combining the optical data 
with the X-ray data does not change the above val- 
ues significantly, though the rising in X-rays between 
To + 50 ks to To + 200 ks is more pronounced than 
at optical/NIR wavelengths. The nearly achromatic be- 
ginning and ending times of the plateau, no significant 
spectral evolution across the break, and nearly same 
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TABLE 3 
Gamma-ray bursts at z > 5. 



GRB 



rodshift 



Reference 



060522 
050814 
060927 
050904 
080913 



5.11±0.01 
5.77±0.12 
5.467 
6.295±0.002 
6.695±0.025 



Ccnko ct al. 



{2006^ 



|Huiz-Velasco et al. N2007 
|Kawai et al~| l |2006| 

this work 



Curran et al. 



GRB 080913 is a proof that GRBs do occur at z > 6.7 
and hence that a mechanism for star formation and evo- 
lut ion leading to a burst of this duration was in place at 
this time. As we move to higher and higher redshifts the 
fraction of massive stars that are of population III (low 
( ^7j7jg|ietallicity) will increase. When this happens is strongly 

dependent on the strength of winds fro m the first stars 

and the e fficiency of mixing of the metals ( Scannapieco et 
al. |2006 |. In most models population III stars and stars 



decay rate before and after the plateau (if averaging 
over the short-scale variability) all suggest that the re- 
brightening is due to the forward shock emission with a 
continuous energy injection, either due to a long term 
central engine injecting energy in the form of Poynt- 
ing flux (Dai & Lu 1998; Zhang & Meszaros 2002) or 
due to spread of the ejecta Lorentz factor distribution 
(Rees & Meszaros 1998; Sari & Meszaros 2000; Zhang 
& Meszaros 2002). In the former model, the Poynting 
flux luminosity from a spinning down central magne- 
tar or black hole is L — Lq(1 + t/To)~ 2 , while in the 
latter model, the mass-Lorentz factor distribution may 
be modeled as M(> T) cx r~ s . The total energy in 
the forward shock during the plateau phase increases 
with time as t a . Using the optical spectral index, we 
can derive a = 2Aa opt /(l + (3 op t) — 1-15 ± 0.16, where 
= 1.22 ±0.15. Therefore, the to- 



*opt 



*opt,l 



*opt,2 



tal energy has increased by a factor of (i&2A&i) a ~ 12 
compared to its initial value Ei. According to Table 
2 of Zhang et al. (2006), we obtain q = 1 — a = 
-0.15 ± 0.16, s = (7a + 3)/(3 - a) = 6.0 ± 0.8 (ISM) 
and s = (3a + 1)/(1 - a) = -29.7 ± 31.8 (wind). The 
value q ~ is quite consistent with the Poynting flux 
injection model for times earlier than the central engine 
spin down time scale (t < Tq). 

3.3. On high-redshift indicators 

Various redshift indicators have been proposed in 
the past with the aim to select, based on high-energy 
data, candidates for high-redshift GRBs. Those include 
Swift /BAT properties like burst duration (Tg > 60 s), 
peak photon flux (< 1 ph/s/cm 2 ), and spectral steepness 
(Campana et al. 2007; Salvaterra et al. 2007; Ukwatta et 
al. 2008), slowly rising GRB profiles (BAT image triggers 
as opposed to rate triggers), or the "pseudo" -redshift (in 
this case z=2. 5-6.0) based on the peak energy (Pelan- 
geon 2006). For GRB 080913, most of these indicators 
have failed (in fact also for GRB 060927 at z=5.47), ex- 
cept for the excess-Nn method (Grupe et al. 2007). As 
these criteria are statistical in nature, outliers are possi- 
ble. However, the parameters of this burst demonstrate 
that one should not bias the search for high-z events on 
those criteria. High-z GRBs are so rare that we cannot 
afford to loose the rare events that do occur. The trigger 
criteria we are setting up for our follow-up observations 
should rather include false triggers than exclude the true 
ones. The only secure way to find these high-z events is 
by simultaneous optical/NIR follow-up observations to 
establish their drop-out nature from the colors of the af- 
terglows. 

3.4. GRBs at z>6 



rmly enriched by population III stars can be found down 
to z > 5. Simulations of supernovae from population III 
stars suggest that the explod ing stars will have hydrogen 
and helium rich outer layers ( Ohkubo et al. 2006 Lawlor 
et al. |2007 1 and hence that they probably will not make 



GRBs at least by the collapsar mechanism. Concerning 
GRBs from binary population III stars it is debated how 
efficient they are in producin g GRBs (Bromm & Loeb 



2006 Belczynski et al. 2007 1 



Detecting high-redshift GRBs with Swift and mea- 
suring their redshifts with ground-based spectroscopy 
is of substantial interest because of the link between 
long-duration GRBs and the star formation rate den- 
sity (SFRD). Indeed, the properties of GRB hosts and 
the distribution of GRB metallicities are consistent with 
the assumption th at GRBs trace the bulk of the star- 
formation at z > 3 



2008a). The existence 



( Jakobsson et al. |2005 Fynbo et al. 
ice of now 5 GRBs at z > 5 (Tab. 



among them two at z > 6 out of a total of ~150 GRBs 
with redshift estimates suggests that the global SFRD at 
z > 5 declines slowly and has a substantial value. Em- 
pirical calibrations of the GRB rate to the SFR at dif- 
ferent redshifts have been attempted, and suggest that 
the high-z SFR must be high to accommodate the sev- 
eral observed high-z GRBs (Chary et al. 2007, Yiiksel et 
al. 2008). Our detection of a GRB at z=6.7 strengthens 
this result. This is also similar to the SFRD in the model 
assumed by Bromm & Loeb (2006) that predicts ~ 10% 
of the Swift GRBs are at z > 5. 

The discovery of bursts like GRB 080913 also allows to 
study, in an unbiased way, some of the main questions in 
the evolution of the Universe: where did most of the star 
formation happen at z > 6, and what was the nature 
of the sources responsible for the re-ionization? There is 
evidence that the bright z > 6 galaxies discovered using 
color-color (drop-out) selection or more advanced pho- 
tometric redshifts are too rare to provide the total star 
formation rate as well as to have done the re-ionization 
(e.g., Bouwens et al. 2007). GRB measurements pro- 
vide the tool to find the more typical galaxies respon- 
sible for the bulk of the p roduction of ionizing photons 
( Ruiz-Velasco et al. ||2007 ), and will allow further study 
of these galaxies in the future (e.g. with JWST or the 
~30 m ground-based telescopes) . 

The re-ionization history of the Universe is currently 
not well constrained by observations. The WMAP data 
(Spergel et al. 2007) and the SDSS quasar data (Fan 
et al. 2006), and the GRB 050904 constraint (Totani et 
al. 2006) can accommodate several distinct re-ionization 
scenarios (e.g. Holder et al. 2003). In order to iden- 
tify the correct scenario, bright beacons in the dark era 
(e.g. z = 7 — 1100) are needed. The discovery of GRB 
080913 above the highest z for QSOs (CFHQS J2329- 
0301 at z = 6.427 ±0.002; Willot et al. 2007) re-enforced 
the possibility to use high-z GRBs to uncover the cosmic 
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re-ionization history. The detection of lensed star forma- 
tion at z > 7 (e.g., Bouwens et al. 2008) suggests that 
star formation in fact took place as early as z ~ 8 — 10 
(at t < 0.63 Gyr). With hints for the first stars having 
formed as early as 20^5 z ^60 (Kogut et al. 2003, Bromm 
& Loeb 2006, Naoz & Bromberg 2007), GRBs are be- 
lieved to exist as early as z ~ 15 — 20, and their gamma- 
ray and IR emissions are bright enough to be detected by 
the current instruments (Lamb & Reichart 2000; Ciardi 
& Loeb 2000; Gou et al. 2004). A caveat has been that 
the GRB host DLAs may have a large column density 
that would bury the signature of IGM absorption, as in 
the case of GRB 050904 (Totani et al. 2006). Cosmo- 
logical SPH simulations (Nagaminc et al. 2008), on the 
other hand, suggest that the GRB host DLA columns 
should degrade at high-z, a prediction dictated by the 
fundamental structure formation theory. The low N(HI) 
associated with GRB 08 09 1 3 ( 5 x 10 20 c m - 2 for the best 
fit value and 1.4 x 10 21 cm~ 2 for the maximum value) is 
generally consistent with such a prediction, though the 
medium N(HI) for z > 4 GRBs (10 2145 cm -2 ) is consis- 
tent with the medium value (10 cm -2 ) for the total 
GRB population (2 < z < 6.3). In any case, high-z 
GRBs offer the promising possibility to probe cosmic re- 
ionization. Brighter GRBs at higher z's with a negligible 
N(HI) would allow mapping the IGM ionization fraction 
as a function of z in the dark era, and hence, greatly 
constrain the possible scenarios of cosmic re-ionization. 

4. CONCLUSIONS 

We have presented the observations that led to the dis- 
covery of GRB 080913 and the recognition of it being at 
rcdshift 6.7, the most distant GRB so far. This discov- 
ery was possible first of all due to the excellent properties 
of the Swift mission with its high sensitivity for GRBs 
and excellent localization capabilities based on the X-ray 
afterglows alone. Furthermore, this study demonstrates 
that 2 m class telescopes are sufficient to properly identify 
and determine reliable photometric redshifts of the most 
distant GRB afterglows. Here it is noteworthy that the 
afterglow of GRB 080913 was several magnitudes fainter 
than that of GRB 050904, the only other z > 6 GRB dis- 
covered so far. Hence, we can probe a significant fraction 



of the luminosity function of afterglows at these redshifts 
with 2 m class telescopes. We also note that many pre- 
viously proposed high-z "indicators" primarily based on 
properties of the prompt emission have failed for GRB 
080913, leaving photo-z determinations the only reliable 
method available. 

While the afterglow of this GRB was particularly faint, 
we note that metallicity measurements as well as the de- 
termination of the neutral hydrogen fraction should be 
possible in general (as for GRB 050904) though likely 
not in all cases. It is important to recognize, however, 
that studying the neutral hydrogen fraction in the high- 
est redshift quasars is hopeless, and for galaxies compli- 
cated by their faintness and m ore structured continuum 
shape ( |McQuinn et al. |2008[ ). Planned instrumentation 
in the near future with higher effi ciency and higher s pec- 
tral resolution such as X-shooter ( Kaper et al. |2008| will 
help substantially in exploiting the promise that GRBs 
have as probes of the re-ionization epoch. 
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Fig. 1. — i! (left), z' (middle) and J (right) band images of the afterglow of GRB 080913 obtained with the 7-channel imager 
GROND at the 2.2m telescope on La Silla / Chile. The circle denotes the Swift /XRT error box. Some of the local standards 
of Tab. [2] are labeled in the i' image; the remaining ones are outside the field shown here. 
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Fig. 2. — Optical/NIR (bottom) and X-ray light curve (top) of the GRB 080913 afterglow. The XMM data point has been converted 
to the Swift/XRT count rate axis using the best-fit spectral model. Several deviations from a canonical power law decay are apparent in 
both, X-rays as well as optical. The early X-ray light curve (To+200 s until To+2000 s) is modeled by two log-Gaussian flares. For the 
optical light curve two different fits are shown: one for the plateau interpretation (dotted lines) and one for the flares interpretation (three 
log-Gaussians; solid line). The first model is motivated by the scenario of a plateau and late energy injection. In the second model, which 
assumes that X-ray flares are associated with optical flares, the last NIR data point is too bright, and consequently needs a third flare. The 
late Gemini upper limit supports this interpretation. The sum of the three log-Gaussian profiles roughly accounts for the flare at X-rays 
at the same time interval. 
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Fig. 3. — Spectral energy distribution of the afterglow of GRB 080913, obtained about 30 min after the GRB. The drop-out shortward 
of the z'-band is clearly visible. This information was used to trigger spectroscopy with VLT/FORS2 which then started with grism 600z, 
selected as compromise between spectral resolution and faintness of the afterglow, about 2 hrs after the GRB. Shifting the later VLT 
magnitudes in BVRIz to the earlier time of the GROND measurement, the upper limits would be B > 22.5, V > 22.6, R > 22.5 and 
/ > 22.0, thus not more constraining than the GROND limits. 
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Fig. 4. — Optical spectrum of the afterglow of GRB 080913 obtained with VLT/FORS2. The spectrum shown corresponds to the 2400 s 
exposure on Sep 13, 2008. The break due to the Lye? forest is clearly visible at 9400 A. The bottom green curve shows sky spectrum, and 
the dotted line superimposed on the spectrum shows the error spectrum (noise) after the sky subtraction. Two model fits are shown (see 
text). 
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Fig. 5. — Combined Swift /XRT and GROND spectral energy distribution for the time interval 200-1800 s after the burst. Shown is the 
fit with a single power law (for parameters see text). 




The optical/NIR afterglow of GRB 080913 (thick black line), shifted to the Rq band and to a redshift 2 = 1, and compared 

At early times, the optical/NIR afterglow of GRB 080913 is fainter 



Fig. 6 

with other GRB afterglows, also at 



than the bulk of the long GRBs with detec 



1 ( Kann et al. 12007b 2008 



cd afterglows, and at late times it is brighter. With the possible exception of GRB 060121 (if it 



lies at z = 4.6, |De Ugarte Postigo et al. |2006[ l, the optical/NIR afterglow of GRB 080913 is also much brighter than even those of bright 
short GRBs. 



